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ABSTRACT 

We present interferometric CO observations of twelve z~2 submillimetre-faint, star- 
forming radio galaxies (SFRGs) which are thought to be ultraluminous infrared galax- 
ies (ULIRGs) possibly dominated by warmer dust (T^ust ^ 40 K) than submillimetre 
galaxies (SMGs) of similar luminosities. Four other CO-observed SFRGs are included 
from the literature, and all observations are taken at the Plateau de Bure Interferom- 
eter (PdBI) in the compact configuration. Ten of the sixteen SFRGs observed in CO 
(63%) are detected at >4<j with a mean inferred molecular gas mass of ~2xl0^^ M©. 
SFRGs trend slightly above the local ULIRG Lfir-L^q relation. Since SFRGs are 
about two times fainter in radio luminosity but exhibit similar CO luminosities to 
SMGs, this suggests SFRGs are slightly more efficient star formers than SMGs at the 
same redshifts. SFRGs also have a narrow mean CO line width, 320±80 kms~^. Many 
SMGs have similarly narrow CO line widths, but very broad features (~900 kms~^) 
are present in a few SMGs and are absent from SFRGs. SFRGs bridge the gap be- 
tween properties of very luminous >5xl0^^ LqSMGs and those of local ULIRGs and 
are consistent with intermediate stage major mergers. We suspect that more moderate- 
luminosity SMGs, not yet surveyed in CO, would show similar molecular gas properties 
to SFRGs. The AGN fraction of SFRGs is consistent with SMGs and is estimated to 
be 0.3±0.1, suggesting that SFRGs are observed near the peak phase of star formation 
activity and not in a later, post-SMG enhanced AGN phase. Excitation analysis of 
one SFRG is consistent with CO excitation observed in SMGs (turning over beyond 
^^C0( J=4^3)). This CO survey of SFRGs serves as a pilot project for the much more 
extensive survey of Herschel and SCUBA- 2 selected sources which only partially over- 
lap with SMGs. Better constraints on the CO properties of a diverse high-z ULIRG 
population are needed from ALMA to determine the evolutionary origin of extreme 
starbursts, and what role ULIRGs serve in catalyzing the formation of massive stellar 
systems in the early Universe. 

Key words: galaxies: evolution — galaxies: high-redshift — galaxies: infrared — 
galaxies: starbursts 
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1 INTRODUCTION 

Ultraluminous infrared galaxies (ULIRGs) exhibit some of 
the most extreme star formation rates in the Universe. The 
volume density of higher redshift ULIRGs peaks at z^2- 
3 (Chapman et al. 2005)— this is also the peak epoch in 



the cosmic star formation rate density and volume den- 



sity of active galactic nuclei (AGN; e.g. Fan et al. 2001 



Richards et al. 2006). Not only does this indicate a pos- 



sible link between supermassive black hole growth and 
rapid star formation, but it also signals the most active 
phase in galaxy evolution and formation. ULIRGs exhibit 
very intense (SFR ^ 2OOM0yr~^), short-lived bursts (r ^ 
100 Myr) of star formation. The possible life cycle of a 
ULIRG, from star- formation dominated, dust-enshrouded 
galaxy, to obscured AGN and then luminous quasar (e.g. 
Sanders et al.|1988 Veilleux et al.|20()9 ), provides a testable 
evolutionary sequence. 

The best studied ULIRGs at high redshift are submil- 
limetre galaxies (SMGs; ,Blain et al.||2002 ) which are char- 
acterised by their detection at 850//m with Ssso ^ 5 mJy. 
While SMGs put powerful constraints on galaxy evolution 
theories and the environments of extreme star formation 
( Greve et al.||2005| [Tacconi et ar||2006| |2008| |Chapman| 



et al. 2005| Pope et al. 2006), their selection is suscepti- 



ble to strong temperature biasing (Bales et al. 112000 Blain 



et al. 2004). At the mean redshift of radio SMGs, z-2.2. 



observations at 850/im sample the Ray leigh- Jeans tail of 
blackbody emission where the observed flux density may 
be approximated by Ssso ocLfir T^^"^^ . Due to the strong 
dependence on dust temperature, the 850/im flux density 
of warm-dust ULIRGs (T^ ^ 40 K) might be much lower 
than cooler dust specimens (T^^ =20-40 K) thus causing the 
warmer-dust galaxies to evade submm detection. This selec- 
tion bias suggests that a large fraction of the 2 ULIRG 
population has not been accounted for in current work on 
high-z star formation. 



Chapman et al. (2004) describe the first observational 



effort to identify warm-dust ULIRGs as a population, 
via the selection of submm-faint radio galaxies (SFRGs) 
with starburst-consistent rest-UV spectra. While they were 
thought to be ULIRGs by their similarities to SMGs (simi- 
lar radio luminosities, optical spectra, stellar masses), with- 
out detection in the far-infrared (FIR), there was no direct 
evidence that their luminosities were in excess of 10^^ L©. 
Ideally, detection at shorter wavelengths in the infrared, at 
A < 500 /xm must be used to confirm a ULIRG 's luminosity 
in the absence of submm detection; Casey et al. ( 2009a ) used 



70 jim. detection to confirm that a subset of the SFRG pop- 
ulation contains a dominating warmer-dust component with 
< Td > = 52 K. While a population of warm-dust ULIRGs 
has been shown to exist, some fundamental questions still 
remain unanswered: are warm-dust ULIRGs in a post-SMG 
AGN heated phase? Could they be triggered by different 
mechanisms than the major mergers said to give rise to cold- 
dust SMGs? 

Investigating the molecular gas content is fundamental 
to the characterisation of star formation properties and gas 
dynamics of a galaxy population. Molecular line transitions 
from carbon monoxide (CO) are a direct probe of the vast 
gas reservoirs that are needed to fuel high star formation 
rates ( Frayer et al. 1999| [Greve et al. 2005 Tacconi et al 



2006| |2008| [Chapman et al.|2 008'). The gas dynamics which 
are derived from these observations shed light on galaxies' 
evolutionary sequences by measuring how disturbed their 
gas reservoirs are and how long they can maintain their star 
formation rates with the observed fuel supply. Recent simu- 
lations work hint that a ULIRG phase may be triggered by 



either major merger interactions (e.g. Narayanan et al. 2009 ) 
or from steady bombardment from low mass fragments (e.g. 
[Pave et al.||201 0); linking observations with these different 
evolutionary scenarios is an essential step in understanding 
galaxy evolution in the early Universe. 

In this paper, we present CO molecular gas observa- 
tions, taken with the IRAM Plateau de Bure Interferometer 
(PdBI), of twelve SFRGs (and four additional SFRGs from 
the literature) to compare the population with SMGs and 
other high redshift star forming galaxies. Section [2] describes 
the sample selection, molecular gas observations and ancil- 
lary data, while section [3] presents our results, in the form 
of derived gas and star formation quantities of the SFRG 
sample. Section [4] discusses the gas properties of the sam- 
ple, compares the population to other high redshift galaxies, 
and hypothesizes on the role of SFRGs in a broader galaxy 
evolution context relative to local ULIRGs and SMGs while 
section [5] concludes. Throughout, we use a A CDM cosmol- 
ogy with Hq 71km s~^ Mpc"\ Qa 0.73 and 0.27 
dHinshaw et al.|2009[ ). 



2 OBSERVATIONS & REDUCTION 

Our sample is drawn from a set of /xJy radio galaxies in 
the GOODS-N, Lockman Hole, Elais-N2, SSA13 and the 



UDS fields using the Chapman et al. (2004) selection of 
submm-faint star forming radio galaxies (SFRGs). They 



were detected in ultra-deep VLA radio maps (Biggs & Ivi- 



2006 



2002 2007 Fomalont et al. 2006) 



Ivison et al 

with Si. 4 GHz ^ 15/iJy at >3cr with an approximate upper 
limit of Si. 4 GHz ^ 1 mJy since strong AGN and radio-bright 
local galaxies were removed from the sample. The //Jy ra- 
dio galaxy population was identified in an effort to isolate 
bright star formers at high redshift, so only the sources with 
non-AGN photometric redshifts of z ^ 1 were included (e.g. 



Chapman et al. 2003). Spectroscopic follow up with Keck 



LRIS revealed starburst spectral features (Chapman et al 



2004 



Reddy et al.[[2006| , mostly at redshifts z > 1. The 



SFRGs with the most reliable spectroscopic redshifts (of- 
ten due to a strong Ly-a emission peak, Gz ^ 0.005) were 
chosen for CO observations at the IRAM Plateau de Bure 
Interferometer. 

We note that all galaxies in our sample satisfy the BzK 
'active' galaxy selection criterion (Daddi et al. 2004), and 
10/14 satisfy the Dust Obscured Galaxy (DOG) selection 
dDey et al.||200 8). All SFRGs have poor rest- UV photome- 
try, so their selection with respect to BX/BM |Steidel et al 
2004) is not constrained. 



Figure [T] shows the distribution in radio luminosity of 
the CO-observed SFRG sample relative to the distributions 
of parent SFRGs, CO-observed SMGs, and parent SMGs. 
We note that the CO observed sample in this paper is about 
two times less radio luminous than the CO-observed SMGs 



which were analyzed in Neri et al. ( 2003 ), Greve et al. ( 2005 ), 



and Tacconi et al. (2006)— an aspect of their selection which 
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Figure 1. The distribution in radio luminosity (and inferred 
FIR luminosity) of CO observed SFRGs (black line-filled) and 
SMGs (gray filled) relative to the distribution in radio lumi- 
nosities of their parent populations: spectroscopically confirmed 
SFRGs (dashed line) and SMGs (dotted hue) without CO obser- 
vations. 



0.9 GHz bandwidth coverage (overlapping polarisations), 
and RGJ123644 which had 1.6 GHz coverage (semi overlap- 
ping polarisations). RGJ123644 had an uncertain spectro- 
scopic redshift and upon initial observations with a 0.9 GHz 
bandwidth, a feature was identified at the edge of the band- 
width and subsequent observations were shifted in redshift, 
thus totaling 1.6 GHz coverage. The total on-source inte- 
gration time varied between ^^4-12 hours. The data were 
processed using the GILDAS packages CLIC and MAP- 
PING and analyzed with our own IDL-based routines. The 
RMS noise of each object's map is also given in Table [l] 
For clarity of presentation, we have re-gridded the data to a 
spectral resolution of 40 MHz for data with >1.5GHz band- 
width and to 20 MHz for data with <1.5GHz bandwidth. 
The bandwidths, binsizes and noise properties of the maps 
are given m Table [1] 

The search for CO line detection was performed by in- 
tegrating over all possible combination of channels at ev- 
ery location in the observed maps using our own iterative 
GILDAS script. From the result we measure the data's sig- 
nal to noise ratio (S/N) at all points in the map, and subse- 
quently investigate signal peaks greater than 4cr. If a detec- 
tion with signal strength >4cr exists within ^" of the target 
position, the observed galaxy is classified as detected. 



traces back to the removal of more luminous spectroscopic 
AGN from the SFRG sample. The equivalent class of spec- 
troscopic AGN are not removed from the SMG sample since 
their detection in the FIR provides sufficient evidence that 
the SMGs are star-formation dominated. This is likely rep- 
resents the most dominant selection bias between the popu- 
lations, however we address more in the discussion in section 
|4.6| We also discuss the AGN fraction of SFRGs at length 
in section [XT] While lower luminosity SMGs have been ob- 
served with PdBI (representing the low luminosity tail on 
the SMGs in Fig[l] Bothwell etal., in preparation), only the 
published CO observed SMGs are included in this paper for 
comparison. 



2.1 PdBI Observations 

Nine SFRGs were observed from June 2008 to October 
2008, while three additional sources were observed in Au- 
gust through October of 2009. The mean redshift of the 
sample is z = 1.9 ± 0.6. 

CO Observations were carried out with PdBI in the 5 
dish D-configuration (i.e. compact). We used the 2mm and 
3 mm receivers tuned to the appropriate frequencies of red- 
shifted CO transitions, as detailed in Table [l] Pointing cen- 
tres were at the VLA positions of the galaxies, and example 
phase calibrators which we used were 0221+067, 1044+719, 
1418+546, 1308+326, and 0954+658 with flux calibrators 
like MWC349, 3C84, and 3C345. The synthesised beam size 
for the configuration varies from ^3-6^^ FWHM. Receiver 
Noise Temperature calibration was obtained every 12min 
using the standard hot /cold-load absorber measurements. 
The antenna gains were found to be consistent with a stan- 
dard range of values from 24-31 JyK~^. We estimate the 
flux density scales to be accurate to about ±15%. 

Data were recorded using both polarisations, offset in 
frequency, covering a 1.8 GHz bandwidth, except in the 
case of RGJ021827, RGJ131207, and RGJ105239 which had 



2.1.1 Offset CO Detections 

Searching for detections within a radius imply that sig- 
nificant positional offsets are acceptable. In our galaxies, 
only three of our sources, RGJ021827, RGJ105209, and 
RGJ105239), have offsets >2". Nominally, an offset >2" 
would be too large to attribute to the targeted source, how- 
ever we note a few factors which can increase the positional 
uncertainty to ^5^^ PdBI D-configuration positional uncer- 
tainty is governed by the source S/N (which is proportional 
to beamsize/2(S/N), ^ 1" ^ but might diverge at small S/N), 
the baseline model uncertainty (^0.25^^ for D-config), and 
the quality of phase calibration (which can translate to 1.5- 
3^^ seeing). Since the S/N of these galaxies is small, then the 
nominal uncertainty likely increases from '-^2^^to . It is 
also possible for the bulk of the gas reservoir in a system 
to be offset from the primary source of starlight, especially 
in the case of major mergers. In addition, we have found 
large positional offsets ^-^5^^ with strong >5cr CO detections 
in the larger samples of SMGs (i.e. CO positional offset from 
radio centroid. Small & Bothwell, private communication). 
Since these detections are centred on the correct redshifts 
and there does not seem to be other radio/UV/IR sources 
corresponding with their positions, we treat RGJ021827, 
RGJ105209, and RGJ105239 as detected; however, we la- 
bel their spectra and CO maps as 'OFFSET' and list them 
separately in Table |2] both as tentatively detected and as 
undetected, giving the upper limits for phase center obser- 
vations. We plot them in subsequent figures as detected, 
although mark them with different symbols to indicate ex- 
clusion from key calculations in our analysis. 



2.1.2 Remarks on Individual Sources 

RGJ163655, RGJ131236, and -^^CO (J=3^2) RGJ123711 
are taken from Chapman et al. ( 2008 ) , the pilot study for the 
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Table 1. PdBI Observation Properties and FIR, Radio Data for SFRGs 
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Table Note s. The top 12 source s were obse rved in our program an d the bottom five CO-observed SFRGs are taken from the 



Daddi et al. 



(20081 and " from Chapman et al. (20081. BW" denotes the bandwidth of observations. Non-detections 



literature: " from 

are 2a upper limits, magnitudes are in AB, and ellipses denote that no data have been taken of that galaxy at the given wavelength 
(24/im, 70/im, 350/im, and 1200/im). RGJ123711 is listed twice, once for its observations taken under our program, and once for the 
observations discussed in Chapman et al. (2008). The 24/im fiux densities of the SSA13 field sources are marked by ^: they are derived 
from IRS spectral observations since no MIPS imaging exists for this field (Casey et al., in preparation). The RMSct is the noise of the 
CO observations averaged over the bandwidth of observations (also over one beamsize at phase centre), while RMSc/i is the noise per 
frequency channel which has width with size Bin, in MHz. The frequency bins are chosen for optimum presentation of the spectra, as 

shown in Figure [3] and are based on bandwidth. The uncertainties on the radio fiux densities are errors on fiux integrated 
measurements and not statistical errors. RGJ021827 was targeted at a redshift of 1.357, which is the redshift of a nearby source; the 
correct UV spectroscopic redshift is 1.361, which is within the bandwidth of our CO observations. RGJ163655 has an Ho; redshift of 
2.192 and a UV redshift of 2.186, both within the bandwidth of the literature observations. The redshift for RGJ131236 has been 
corrected from earlier measurement of 2.240. Its CO observations were taken at 2.240, so unfortunately the CO[3-2] line at 2.224 falls at 
the edge of the PdBI bandwidth, which makes it difiicult to put constraining limits on its CO emission. 



sample observed in this paper. These data, when presented 
in this paper, are different from the spectra and maps pre- 
sented in Chapman et al. due to a different reduction. We 
decided to re-reduce their data to be consistent with our 
reduction, making significant improvements in both phase 
and amplitude calibration. The differences do not change 
the results of Chapman etal., but a few changes in minor 
conclusions are noted later on in the results and discussion. 
The Chapman etal. CO observations of RGJ131236 were 
carried out assuming z =2.240 however more recent rest- 
UV spectroscopic observations indicate a different rest-UV 
redshift, z =2.224. The final channel in the data cube might 
suggest a flux excess at the correct redshift, although we 
cannot distinguish between a bright line cutoff at the edge, 
a broad faint line, continuum, or noise spike. We treat this 
source as undetected with the associated noise properties 
of our observations. We also note that RGJ123711 was ob- 



served as part of [Chapman et al.| ( |2008 ) in ^^CO(J=3^2), 
but newer, higher S/N observations in ^CO(J=4^3) were 
taken as part of this program. For comparison, we include 
the Chapman et al. results in Tables [l] and [2] In the re- 
reduction of RG J 163655 observations, we find a possible 
companion CO source at the same redshift offset by ~8^^ 
to the northwest which is not seen in the Chapman et al. 



results, however significant improvements on phase and am- 
plitude calibration have been made since. This source is 
discussed briefly in section |3.1| The new data taken for 
RGJ123711 reveals a marginaP^CO( J=4^3) detection for 
the SMG SMMJ123711. 98+621325. 7 (also called HDF255); 
its integrated CO flux (detected at 4.5cr) is /i2co(j=^->^3) — 
0.54ib0.12 Jy km s~^ ; its spectrum and properties will be dis- 
cussed more at length in a paper summarising CO properties 
of observed SMGs to date (Bothwell etal. in preparation). 



2.2 Archival Observations 

This paper also includes the CO-observations of two addi- 



tional SFRGs from the literature. Daddi et al. 



(2008) an- 



alyzed RGJ123710 and RGJ123626 as being high redshift 
normal spiral galaxies which have very low star formation 
efficiencies (a study expanded upon in Daddi et al. 2010). 



The active BzK galaxies which Daddi etal. survey have 
significant overlap with the SFRG population, since their 
CO sources require radio detection, thus they have ULIRG 
luminosities implied from the radio. Both RGJ123710 and 
RGJ123626 are selected as SFRGs via the [Chapman et al.| 
( 2004| ) method and are likely very-luminous star formers, 
with higher SFRs than most active BzKs. RGJ123710 has 
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also been detected at TO/xm (Casey et al. 2009a), directly 



confirming that it is a ULIRG with a warm dust temper- 
ature. Interpreting the BzK active galaxy population as 
ULIRGs rather than "high-z normal spirals" is perhaps sen- 
sible for these reasons. 



2.3 Multiwavelength Data 

Radio fluxes from VLA B-array maps (5^^ resolution) are 
taken from Richards' (»2000) and fMorrison et al.l (|2008|), 



where the uncertainty represents the error in the ex- 
tracted flux measurement rather than statistical error. 
High-resolution observations from the Multi-Element Ra- 
dio Linked Interferometer Network (MERLIN; |Thomasson| 
1986) were obtained for the sources in GOODS-N as de- 
scribed in Muxlow et al. ( [2005J , and a combined MER- 
LIN+VLA map was constructed with an RMS noise of 
4.0 /iJy beam~^. A similar map was constructed i n Lockman 
Hole with an RMS noise of 4.5 /xJy beam~^ (Biggs & Ivisonl 



2008). The combined MERLIN+VLA maps have positional 



accuracies of tens of mas and restoring circular beam sizes 
of 0.4'' (GOODS-N) and 0.5''(Lockman Hole). 

RGJ105209, RGJ123711 and RGJ131208 were all ob- 
served at 350/im with the SHARC2 camera (Dowell 
ah] [2003 ) on the Caltech Submillimeter Observatory 



et 



and reduced with CRUSH software ( |Kovacs et al. 2006). 
With on-source integration times of 950s, 1080s and 200s, 
RGJ105209, RGJ123711 and RGJ131208 were have flux den- 
sities of 19±7mJy, 6ib9mJy, and 48ibl7mJy respectively. A 
fourth SFRG, RGJ163655, was observed at 350/im with the 
CSO previously ( Chapman et al.||2008 ). None are detected 
at >3cr. While the acquisition of ground-based 350 /xm ob- 
servations is laborious and dependent on the best weather 
conditions (tcso < 0.06, measured at 225 GHz), the relative 
depth of our observations is comparable to expected confu- 
sion limits for the Herschel Space Observatory at 350 /xm. 
While Herschel will detect >40 mJy 350 /im sources at high- 
z regularly, sources with flux density ^ 20 mJy will need 
follow up from ground-based facilities like the CSO for more 
precise flux density estimates and better SED constraints. 

The 1200/im flux limits for GOODS-N and Lockman 
Hole come from the Max-Planck Millimeter Bolometer 
(MA MBO, with a mean RMS of ^0.8 mJy; |Greve et ^1.] 
' 2008| ) and 850/im flux limits from the Submm Common User 
Bolometric Array (SCUBA, with a mean RMS of ^1.6 mJy; 
Borys et al.|2003| [Coppin et al.|2006t . All fields are covered 
with Spitzer IRAC (3.6, 4.5, 5.8, and 8.0/im) and MIPS 
(24 and 70/im), however at greatest depths in all bands 
in GOODS-N. Optical photometry in GOODS-N is from 
the HST AC^ using the F435W, F606W, F814W, and 
F850LP filters (B, V, i and z bands). The Lockman Hole 
has HST ACS F814W (PI: Chapman HST 7057) in addition 
to extensive optical photometry from Subaru/Suprime-Cam 
( Miyazaki et al.||2002 ). X-ray fluxes are measured from the 



^ Based on observations made with the NASA/ESA Hubble 
Space Telescope, and obtained from the Hubble Legacy Archive, 
which is a collaboration between the Space Telescope Science In- 
stitute (STScI/NASA), the Space Telescope European Coordinat- 
ing Facility (STECF/ESA) and the Canadian Astronomy Data 
Centre (CADC/NRC/CSA). 



Chandra/XMM maps of GOODS-N ([Alexander et al.|2003|), 
Lockman H ole ( [Brunner et al.||2008| ), SSA13 ( [Mushotzky 
eraL][2Q00t , and SXDF ( jUeda et al.f2008 ). 



2.4 RGJ105209 

Recent work on mid-infrared spectra of SFRGs from the 
Spitzer InfvaRed Spectrograph (IRS) (Casey et al., in prepa- 
ration) has revealed that the 24/im source that lies at the 
position of RG J 105209 (at the radio position) has a domi- 
nant PAH redshift of 2.37, in contrast to the CO and UV 
spectroscopic redshift of 2.112. We suspect that two differ- 
ent systems overlap to create these discrepant redshifts, but 
note that only one source is visible in the rest-UV. The high- 
resolution MERLIN+VLA radio imaging also shows only 
one point source slightly offset 1.2^^ from the UV source. 
Due to positional uncertainties, it is difficult to pinpoint 
which source is generating the radio emission. For that rea- 
son, we assume in this paper that the radio is associated with 
the 2.112 CO source. We note that similar positional overlap 
phenomena have occurred with star forming galaxies before, 
for example there are several examples of sources with two 
distinct redshifts for one continuum source in the surveys of 
Steidel et a L](r2004) andl Reddy et al.| ( |2006l |2008| ), and the 
density of bright star forming z~2 sources is high enough 
that overlap will occur with non-negligible probability. 



3 ANALYSIS AND RESULTS 

3.1 Derived Properties from CO Data 

Table |2] lists the integrated line fluxes and limits for all 
SFRGs in our sample. Out of the twelve sources in our ob- 
serving program, seven SFRGs were detected at ^ 4a. In- 
cluding the literature SFRGs, 10/16 are CO detected, and 
one 'undetected source' has insufficient data to determine 
detection (RGJ131236, with an ambiguous redshift). Three 
of our detected sources have significant positional offsets 
from their radio positions, thus we have classified these as 
'tentative' detections in Tableland excluded them from cal- 
culations in our analysis. The detection fraction of SFRGs, 
10/16 (63%, including offset sources), is similar to the de- 



tection fraction of SMGs in Greve et al. ( 2005 ) and Coppin 
et al.| ( |2008] ). 

While we do not expect continuum fiux from any 
SFRGs, we test for it based on FIR detection limits. To 
test the detectability of blackbody continuum radiation, we 



interpolate the best fit FIR SEDs (see section 3.6) to esti- 
mate the fiux density at the frequency of PdBI observations 
(zy^bs = 80 - 150 GHz). An SFRGs in the sample have esti- 
mated continuum fiux densities <0.13mJy, which is of order 
the RMS noise in our CO spectra, thus not detectable. 

Each remaining SFRG spectrum was fit with a three 
parameter Gaussian using a least squares fitting algorithm, 
where the spectral noise is simulated as a function of the 
RMS channel noise (see [Coppin et al.|[2007 for details). 
Since no detectable continuum emission is expected we do 
not fit line profiles with constant baselines. Figure [2] shows 
the velocity- averaged spatial maps of the mm line observa- 
tions (with VLA centroids marked with crosses) and Fig- 
ure [3] shows the extracted spectra of the CO detections with 
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Figure 2. CO maps of SFRG sample in AO" x 40" cutouts, sorted first by detected CO (the first eight are detected), then by RA. 
For those galaxies detected in CO, the maps are integrated over the optimal velocity channels corresponding to the detected CO line. 
Black contours are integer multiples of the RMS starting at 2a and thin gray contours are negative integer multiples of RMS starting at 
-2(7. The relative fiux scales on the maps are indicated by the color-bars to the right. The spectra in Figure [3] are extracted within one 
beamsize (white contours in the lower right, outer contour is FWHM) centred on the highest S/N point in the map (large cross). The 
three sources at top are labelled "Offset" since their CO peak is ^3-5" offset from phase center, a greater offset than would be caused 
by positional inaccuracies in the instrument. For galaxies not detected in CO, the maps are integrated ^jver ^^J^^^^^j^J^^^^qqq^^ 
spectra are extracted from the map center (small cross) . W 5 ? |_ 



CO in z>l submm-faint ULIRGs 



REDSHIFT OF C0[2-ll 
1.350 1.355 1.360 



REDSHIFT OF CO[3-2l 
2.10 2.11 2.12 



REDSHIFT OF C0[2-l] 
1.820 1.825 




; RGJ 105209 




OFFSET ; 








-1500 -1000 -500 
VELOCITY (km s"') 

REDSHIFT OF C0[3-2] 
2.08 2.09 2.10 



-1000 1000 

VELOCITY (km s"') 

REDSHIFT OF C0[2-l] 
1.42 1.43 ^ " 



L RGJ 125644 * 










-1500 -1000 -500 500 1000 1500 
VELOCITY (km s"^ 

REDSHIFT OF C0[4-3] 
1.980 1.985 1.990 1.995 2.000 2.005 2.010 



A 1 




-1000 1000 

VELOCITY (km s"') 



-1000 1000 
VELOCITY (km s~') 



-1500 -1000 -500 500 1000 1500 
VELOCITY (km s"') 



REDSHIFT OF C0[3-2] 
2.22 2.23 2.24 2.25 



REDSHIFT OF CO[3-2] 
2.180 2.185 2.190 2.195 2.200 2.205 2.210 



2.5 F 
2.0 
1.5 
1.0 



0.0 
-0.5 



RGJ151208 



h 



t 



2.0 
1.5 
1.0 
0.5 
0.0 
-0.5 
-1.0 



RGJ 16. 

m 


5655 


y 




i 


1 — 






JT 


r 



-1000 1000 

VELOCITY (km s"') 



-1500 -1000 -500 500 
VELOCITY (km s~') 



REDSHIFT OF C0[4-3] 
4 3.66 3.B8 



REDSHIFT OF C0[2-1] 
1.27 1.28 



RGJ1 23642 




It 




A 




REDSHIFT OF C0[2-1] 
1.50 1.51 1.52 



RGJ125718 



n 



1^ 



-2000 -1000 1000 
VELOCITY (km s"') 

REDSHIFT OF C0[2-ll 
1.47 1.48 1.49 1.50 



-1000 1000 
VELOCITY (km s"') 

REDSHIFT OF C0j"2-l] 
5 1.530 1.535 



-2000 -1000 1000 2000 
VELOCITY (km s"') 

REDSHIFT OF C0[3-2] 
2.230 2.235 2.240 2.245 2.250 



2.0 
1.5 
1.0 
0.5 
0.0 
-0.5 
-1.0 



; RGJ123707 






















RGJ131236 




-2000 -1000 1000 
VELOCITY (km s~') 



-500 500 
VELOCITY (km s~') 



-500 500 
VELOCITY (km s~') 



Figure 3. The CO spectra for SFRGs. Spectra are extracted within one beamsize centred at the highest S/N peak for detected sources 
and at phase center for non-detections (i.e. the galaxies' radio positions). Vertical dashed lines indicate the galaxies optically derived 
redshifts. Gray brackets highlight the channel range integrated over to produce the maps in Fig[2] The first three sources labeled "Offset" 
are detected at large offset positions ~3-5^^ from phase center. At bottom right, we show the spectrum for RGJ131236, the SFRG observed 
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Table 2. Derived Gas Properties of submm- faint ULIRGs 
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Table Notes. The observed CO properties of the sample. The top seven sources have clean >4cr detections of CO at phase center, the 
next three sources (labeled 'Offset') have CO detections at the correct redshift but offset ~3-5^^ from phase center, and the remaining 
six sources are undetected in C O. -^^p and \co cire given for the observed CO transition, which does not depend on source excitation. 
Then assuming the Weiss et al. (2007) SMG excitation ladder (see text for description), we convert to /co[i-o] ^i^d ^'QQyi_Qy which is 
given in the subsequent columns. A\/co is the FWHM of the fitted feature shown in Fig. [2] The SFRG with a double peaked feature 
(RGJ123711) has the details of each feature listed separately. Ico 2-cr limits for undetected SFRGs are calculated assuming a 
AVco =320 kms-i (the mean FWHM of the detected-SFRG sample). 



overlaid Gaussian fits and the original optical spectroscopic 
redshifts. The 2D maps for CO-detected SFRGs are inte- 
grated over the channel range which produces the highest 
signal to noise line profile in the spectrum. The 2D maps of 
the undetected sample are integrated over all velocity chan- 
nels. CO detected SFRGs have CO redshifts which agree 
with their optical spectroscopic features. Every spectrum 
was also tested against two three-parameter Gaussian fits 
(i.e. a double peaked line) in a x squared goodness of fit 
algorithm. Only one galaxy's line fits were significantly im- 
proved by fitting to a double peaked Gaussian. The detection 
of RGJ 123711 is double peaked with one component centred 
at zco = 1.988 and another at zco = 1.996 which bracket its 
optical redshift (z = 1.995). The line properties of both com- 
ponents of RGJ123711 are listed separately in Table [2] the 
total system is represented by the sum of the two separate 
gas masses and sum of dynamical masses (we treat them 
as separate systems since the peaks are asymmetrical and 
more likely indicative of a merger than rotating disk). The 
ratio for double-peaked CO in SFRGs is thus 1/8?^ 13%, 
comparable to 5/30^^17% for SMGs. The weighted mean 
line width of this sample, excluding 'offset' detections is 



320±80 km s~^ (note that it increases to 370±110 kms~4f 
the offset sources are detections). 

The 2(7 line intensity limits for the CO-undetected 
SFRGs are calculated using the following: 



/co < 2RMSch (AVco^i^) 



1/2 



(1) 



where RMSch is the channel noise from Table ^ and dv is 
the spectral resolution, the bins ize value given in Table [l] 
converted to kms~ fas in Boselli et al. 12002). These limits 



are given in Table 
SFRGs (in Figure 



2] The spatial maps for the undetected 
2) are integrated over the entire band- 
width of observations. None of the maps are corrected for 
primary beam attenuation; given the 40x40^^ size, the PBA 
correction is <2 even for all peripheral areas of the map. 

The spatial map for RGJ163655 shows another bright 
companion CO emitter ~8'^ to the northwest of RGJ 163655 
not previously identified in the reduction of Chapman et al. 
2008. Extraction of its spectrum reveals a line centred on the 
same frequency as the RGJ163655 ^^C0( J=3^2) Une, how- 
ever there are no radio sources or spectroscopically identified 
sources nearby. If followup CO observations of RGJ 163655 
reveal the same feature in different CO transitions, then 
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high resolution CO will be needed to precisely measure the 
companion galaxy's position, and thus its multi-wavelength 
properties. 

A lack of multiple CO transition data has meant that 
many past studies have relied on the assumption of constant 
brightness temperature to convert between higher transi- 
tion line luminosities and Lco(i-o)- However, recent obser- 
vations of multiple CO transitions in high redshift galaxies 
has shown that the transitions are often not in thermal equi- 
li brium (.Dannerbauer et al.|2009t . The three z ^ 2.5 SMGs 
of Weiss et al. ( 2007 ) are consistent with being in thermal 



equilibrium up to a J=3 transition yet not towards higher- 
J values, while there are several other high redshift sources 
which are much more excited, for example, FSC 10214+4724 



( Scoville et aL]|l995| ) and APM 08279+5255 ( |WeiB et aL 
2007) which have nearly constant brightness temperatures 
out to J=6) or much less excited (ERO J16450+4626 
turns over at CO (3-2)). In this paper we adopt the bright- 
ness temperature conversions inferred from the spectral 
line energy distributions for three SMGs observed in |Weiss| 
et al. ( 2007| ; explicitly, we assume the following line flux 
ratios derived from LVG models: Sco2i/Scoio = 3.0+0.1, 
Sco32/Scoio = 6.8+0.5 and Sco43/Scoio = 10.0+0.8. One 
SFRG has multiple line data (RGJ123711) and can be an- 
alyzed for its excitation. This work is presented in section 

E2 

More recent work (Ivison et al. 2010c Danielson 
et al.||2010| ) on ^^CO(J=1^0) observations in SMGs show 
that the assumption of a constant brightness tempera- 
ture can lead to underestimations of gas mass by factors 
of ^2x. They find that the ratio of luminosities between 
i2cO(J=3^2) and ^^CO(J=1^0) is r3-2/i-o =0.55 and 
^^0.67 respectively. Our excitation assumption implies a ra- 
tio r3_2/i-o ^^0.78, slightly higher than these recent stud- 
ies, but well within the uncertainties of the *S'i2co(j=3->2) 
measurements of Ivison etal. and Danielson etal.. The de- 
rived line luminosities for our observations are listed in Ta- 
ble |2] in their observed transition as well as in our conversion 



to Lr 



CO(l-O)- 



3.2 CO Excitation of RGJ123711 

The observations of RGJ123711 in ^^00(7=3^2) and 
^^CO(J=4^3) allow an analysis of the CO excitation in 
the system. Figure [4] shows the spectrum of both lines and 
the proposed CO spectral line energy distribution (SLED) 
for the source; since the ^^CO(J=3^2) observations do not 
cover the entire velocity range of the emitting gas, a direct 
comparison of the line strengths can only be done in the 
channels where data overlap (between -700 and 300kms~^) 
marked by vertical dashed lines in Figure |4] The ratio of 
fiuxes is then Sco[4-3]/Sco[3-2] = l-6±0.4. We use Monte 
Carlo testing to estimate the added uncertainty introduced 
by not including the missing portion of ^^C0( J=3^2) (from 
-1150 to -700 kms""*^) and we remeasure this ratio to be 
Sco[4-3]/Sco[3-2] = 2.0+0.8. This agrees with the mean [4- 
3]/[3-2] ratio observed for SMGs, 1.4+0.2 ( [Weiss et al.|2007t 
and the expected ratio if constant brightness temperature is 
assumed, 1.8, however it appears to be inconsistent with the 
ratio for the Milky Way, -0.7 ([Fixsen et aT]|1999). 



We contrast this result with those of Weiss et al. (2007), 
Dannerbauer et a L|( |2009| ), and |Danielson et aL "(2010) in the 



right panel of Fig. [4] Although we have taken measurements 
for J =3 and 4 only and other work is restricted to J <3 
transitions, we see that RGJ123626 and RGJ123711, the 
two SFRGs included here, seem to exhibit different excita- 
tion levels where RGJ 123626 is more consistent with Milky 
Way type excitation and RGJ 1237 11 is consistent with the 
Weiss et al.| ( |2007[ ) SMGs. The issue of gas excitation is 



essential to the physical interpretation of high-z star form- 
ers, and more/multiple line transitions are necessary for a 



CO transitions, like ^^CO(J=1^0) in 


Riechers et al.||2006 


Mainline et al.||2006 Ivison et al.|2010c 


Carilh et al.|2010). 



3.3 Gas Mass and Dynamical Mass 

To derive H2 masses we assume the ULIRG H2 / CO gas con- 
versi on factor X Mh^/L'cq 0.8 M© (Kkms"Vc^)"^ 
from | Downes &: Solomon| ( |1998 ). It is worth noting, how- 
ever, that the CO to H2 conversion factor varies substan- 
tially depending on the type of source being sampled, and 
it can range from 0.8 in ULIRGs to 4-5 in normal spi- 
ral, late-type galaxies (see discussion in the Appendix of 
The basis of our choice oi X — 
0.8 M0 (K kms~ pc )~^ is based on the merger, or dis- 
turbed nature of ULIRGs and their high star formation 
rates. Mergers are confined by ram pressure which pervades 
the whole system, while quiescent discs are confined gravi- 
tationally leading to gas fragmenting into clouds (per the 
jeans mass) and thus much higher ratios of H2 mass to 
unit CO. We note that higher gas conversion factors are 
possible and perhaps likely if SFRGs were to represent an 
intermediate stage or 'less extreme' population of galax- 
ies than SMGs. For this reason, we advise that our H2- 
dependent quantities (like Mh2, gas fraction, etc) be taken 
with a hint of caution, since the full range of galaxy prop- 
erties and conversion factors (from Xqo =0.8-4.5) sug- 
gest variations of Mn^ of order 0.75 dex and variations of 
gas fractions of order 0.6. We also include gas mass esti- 
mates for Xqo = 4.5 M0 (K kms""*^ pc^)""*^ for contrast in 
Table |3] but proceed with Xco=0.8 for our analysis. Us- 
ing Xco=0.8 for ULIRGs, the mean gas mass of the sam- 
ple is (2.1+0.7) X 10^° Mo, which is roughly h alf the mean 
H2 mass of CO-observed SMGs (5.1x10^° M© 
2003J 



Greve et al. 



2005 Neri et al. 



Dynamical mass is dependent on the galaxy's incli- 
nation angle; we use the average inclination correction of 
(sini) =1/2 (in other words i = 30^, corresponding to a 
random distribution in galaxy angles between and 90°). 
Thus 



dyn SZTl % 



Ca^r 
G 



4.215 X l(fC/Wcor 



(2) 



where C —A (the value adopted for mergers; Genzel et al 
2003) and r is given in kpc and AVco in kms~ . We note 
however that the difference between assuming a merger {C — 
4), spheroid (C = 5), or rotating disk (C = 3.4) is on the 
order of the uncertainty in MdynS^^^ i caused by uncertainty 
in line width. 

Deriving a dynamical mass is dependent on a resolved 
spatial size measurement (r in Eq 2) , which is unfortunately 
unavailable with this low resolution, 3-6^^ beam size (D- 
configuration PdBI) data. In place of measuring CO sizes 
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Figure 4. Left: Both -'^^CO( J=3— ^2) and ■'^^CO( J=4— ^3) observed transitions for RGJ123711. Determined from the spectral hnes seen 
in the higher S/N CO [4-3] observation, the velocity range of the emitting region which overlaps in each data set is marked by two vertical 



dashed lines at -700 kms and 300 kms . Right: the CO SLED normalised to I 



CO[3-2] 



for RG J 123711 overlaid with literature results 



from other high-redshift sources (Weiss et al. 2 007||Danielson et al. 2010) and the Milky Way ( jPixsen et al . 1999). One of these sources 
is RGJ123710, an SFRG observed by |Dannerbauer et al.| (|2009|>. We also include LVG model no. 1 from Dannerbauer etal., representing 
a low-excitation Milky Way type source. By normalising to Ico[3— 2]5 we illustrate that it is difficult to draw conclusions on the nature 
of the CO SLED without observing at least three CO transitions, particularly in the region where the SLEDs are suspected of turning 
over near CO [4-3]. 



explicitly, we assume that the SFRGs of this sample have 



gas emitting regions roughly the size of SMGs ( Tacconi et al. 
2008| which have R1/2 =2zblkpc. This is supported by our 



analysis of high-resolution MERLIN+VLA radio imaging 



discussed in the next section (in section 3.4), where SFRGs 



and SMGs are shown to have similar radio sizes and SFRG 



sizes average to 2.3=1=0.8 kpc. Bothwell et al. (2010) have 



taken high-resolution CO observations of RGJ123711 and 
found its size to be consistent with this assumption (its CO 
emission extends over 12.6kpc in one direction and ^5kpc 
in the other, corresponding to a 3.9kpc effective radius, only 
0.2kpc larger than our measured MERLIN+VLA radio size). 
RGJ123711 is perceived to be a merging system that has two 
components of roughly equal spatial extent ^ 2kpc, thus its 
effective size is about ^2x larger than the rest of the galaxies 
in our sample, although we still use the MERLIN+VLA sizes 
as priors for CO size due to lack of resolved CO data. For 
the sources with measured MERLIN+VLA sizes, we use the 
i^eff as the radius r in Equation 2, and for sources without 
MERLIN+VLA data, we assume that r is equal to the mean 
of the measured MERLIN+VLA sizes, r = 2.3+0.8 kpc. 

It is important to note that the extent and morphol- 
ogy of the MERLIN+VLA radio emission can affect the in- 
terpretation of the dynamical mass estimate. Three objects 
have large radio/UV offsets or very extended radio emis- 
sion: RG J 105209 (whose offset is discussed in section 2.4), 
RG J 123707 (extended across a ^7 kpc x 17 kpc area), and 
RGJ123710 (two knots separated by --8kpc). If the CO gas 
were to trace the MERLIN+VLA morphology perfectly, as is 
our a priori assumption, then the radius of each knot of radio 
emission would be used to estimate dynamical masses and 
they would be summed. Depending on the distribution of 



these knots, the true dynamical mass would be calculated by 
considering the size of each knot and the distance separating 
them, however, within the uncertainty of our measurements 
and assumptions (e.g. inclination angle, line-width uncer- 
tainty and uncertainty on the value of C) we consider our 
approach of circularising the MERLIN+VLA sizes and infer- 
ring a radius as accurate. We caution that dynamical masses 
are potentially slightly underestimated, due to the size and 
distribution of the CO gas relative to radio emission. For 
example, radio sizes are potentially inconsistent with CO 
sizes, as with the B-configuration data of RGJ123626 and 



RGJ123710 presented in Da ddi et aL] ( |2010t , which suggest 
larger dynamical mass estimates by a factor proportional to 
the radius. 

When applied to the whole SFRG sample, 
we derive a median SFRG dynamical mass of 
Mdyn =(7.2l^:^)xl0^° MQ/{sin^i) 2.9x10^^ M©. 

The same value of C is used to compare 
with mean dynamical mass of SMGs, which 
is Mdyn(6'MG) =(1.5lJ:^)xlO^^ M0/(sm2z) ^ 
6.0x10^^ M©. 



3.4 MERLIN+VLA Radio Morphology 

We use the high-resolution MERLIN+VLA radio maps to 
map the extension of starburst emission and assess the con- 
tribution of AGN by considering their radio morphology, 
which is shown as contour overlays on optical imaging in 
Fig. [5] With resolutions of 0.4^^ and 0.5^^ per beam (for 
GOODS-N and the Lockman Hole respectively), the small- 
est resolvable structure at z~1.5 would be 4-5 kpc across. 
The typical size of an AGN emission region in a /iJy radio 
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Figure 5. Contours of MERLIN+VLA high-resolution 0.4-0.5^^ radio emission overlayed on ?)" x ?)" optical imaging for the CO-observed 
SFRGs. The optical imaging of GOODS-N sources is from H ST- ACS B, V, and z bands (presented in tricolor). The Lockman Hole 
sources have ACS i-band (black and white cutouts). The MERLIN+VLA beam shape is shown in the lower right of each panel in white 
(outer contour is the FWHM). The contour levels plotted integer multiples of a starting at 3cr (except for RGJ123711 and RGJ123653, 
where only 3, 5, and 7a contours are shown for clarity). Structure in radio emission is seen in these galaxies on 1-8 kpc scales. 
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galaxy at z > 1 (see [Casey et al. 2009b ) at radio wavelengths 
is much less than Ikpc, implying that an AGN dominated 
source would be completely unresolved in MERLIN+VLA 
radio maps (e.g. as seen in the high stellar mass, giant el- 
liptical systems of Casey et al.|[2009b ). Fig. [s] shows that 
each of these galaxies has extended emission regions on large 
scales irregularly spread across large regions of the galaxy, 
suggestive of spatially distributed star formation, unlikely 
to be generated by AGN. About 10-60% of the total flux 
from the source is estimated to be resolved out by the high- 
resolution imaging, which is dependent on the galaxies' ex- 
tended or compact morphologies. We measure an effective 
radius of the star forming area by isolating the regions where 
MERLIN+VLA radio emission is significant to >3cr, and we 
then take the square root of the surface area over pi (effec- 
tively circularised), Reff, which is given in Table |3|. While 
the morphologies are irregular (and are not in fact circular) , 
we find that the effective radii average to 2.3±0.8kpc, which 
agrees with the size measurements of SMGs in [Biggs Ivi-| 
' ' ([2008t and [Chapman et al.j (|2004|. We use the agree- 



The FIR luminosity is then converted into a star for- 
mation rate (SFR) using the following relation: 



ment of MERLIN+VLA SFRG and SMG sizes to partly 
justify our assumption of similar CO sizes between the pop- 
ulations. 

Since as much as 60% of the radio flux can be resolved 
out in the high-resolution maps, we caution that this effec- 
tive radius might be an underestimate, and that it could 
increase by factors of 1.1-1.5 x if all of the flux is accounted 
for; however, this adjustment is not made to our effective 
radii measurements since it is highly uncertain and relies on 
assumptions of the distribution of star formation activity in 
the outskirts of each source. Similarly, this underestimate 
of effective radius would propagate to the derivation of dy- 
namical masses. 



3.5 Star Formation Rates, Densities and 
Efficiencies 

We estimate SFRs from VLA radio luminosities, using the 
radio/FIR correlation for star forming galaxies (e.g 
et al.[[1985[ |Condon[[T992l [Sanders Mirabel|1996t : 

(l-a) 



Helou 



Lfir = (3.583 X 10"^^) Di Sia (1 + 



(3) 



where is luminosity distance in cm, Lfir is evaluated 
from 8-1000/im and is in L©, 5'i.4 is the radio flux density 
at 1.4 GHz in units of /xJy. A Salpeter initial mass function 
is assumed. The factor of (l + z)^^"""-* accounts for band- 
width compression and the radio K-correction to the rest- 
frame luminosity; a is the synchrotron slope index, here 
taken to be 0.75 (|Yun et al. 2001). Since this is an em- 
pirical relation defined by the mean ratio q between FIR 
and radio luminosities, there are a few different versions in 
use in the literature. This relation uses the bolometric ratio 
between FIR flux and radio flux of giR = 2.46 (e.g. Ivison 
et al.|[2010a|b Casey et alT[|2010 ). The relation is tradition- 



ally defined by g = 2.35 ([Sanders Mirabel|p^96 ) and was 



used by [Chapman et al. [( [20Q5f to calculate FIR luminosities 
for the redshift-identified SMGs (only differing by factors of 
20% between different assumed q ratios). In most cases, the 
FIR luminosities derived from the radio are consistent with 
the submm detection limits and any flux density measure- 
ments at shorter wavelengths in the FIR (e.g. at 70/im or at 
350/im). 



SFR{MQyr-^) = 1.7 x 10"'° Lfir ( L^ 



(4) 



from Kennicutt (1998). The inner quartile (25-75 percent) 



of star formation rates is 4OO-16OOM0 yr~ with median 
7OOM0yr~^. Both derived quantities, Lfir and SFR, are 
given m Table [3] We also derive UV-inferred SFRs usmg 
v-hand and i-band optical photometry. Since all SFRGs are 
quite optically faint {i > 22) and FIR luminous, they are po- 
tentially subject to significant (yet uncertain) extinction fac- 
tors. While not taken into account for the UV-inferred SFRs 
themselves, the large extinction factors may be deduced 
from the large disparities between SFRradio and SFRuv 
given in tablejsjwith a typical ratio of SFRradio /SFRuv = 30 
(we do not consider SFRradio /SFRuv to have AGN contam- 
ination based on the previous section). 

Additional evidence that extinction caused by dust is 
significant comes in the comparison between rest-UV mor- 
phologies and their MERLIN+VLA radio morphologies as 
seen in Fig.js] Since the brightest and/or bluest components 
of the rest-UV emission and radio emission do not often 
coincide, it indicates that the dustiest, most FIR-luminous 
regions in the galaxy could be highly obscured. 

Total star formation rate densities, Ssfr, are estimated 
by dividing these SFRs (from VLA fluxes) by the surface ar- 
eas of MERLIN+VLA radio emission regions. We measure 
an effective star formation surface area by isolating the re- 
gions where MERLIN+VLA radio emission is significant to 
>3cr (the values for effective radius of radio emission and 
SFR density, Ssfr, are given in Table |3|. The median to- 
tal SFR density is 30^20 yr"^ kpc"^ Note that this is 
potentially an overestimation of the SFR density since the 
radio flux contained within the >3a MERLIN+VLA re- 
gion constitutes a fraction of the total VLA radio flux (it 
is inferred that this fraction is within 40-90%). For this rea- 
son we also compute alternate star formation rate densities 
using the total integrated MERLIN+VLA flux within the 
>3(7 emission area. We convert the flux to a radio luminos- 
ity, then FIR luminosity, then SFR, and divide it by the 
area. The median value for this alternate SFR density is 
lOl^Q M© yr~^ kpc~^ which is on average, three times less 
than the total Ssfr- 

Comparing either SFR density measurements to their 
theoretical maximum— the maximum gas density divided by 
the local dynamical time (see equation 5 of Elmegreen|1999 



we use tcZyn=4x 10 Vi")— we can determine if the implied 
SFR density exceeds the theoretical prediction. While local 
ULIRGs with Esfr ~200 M© yr~^ kpc~^ are forming stars 
at their theoretical maximum (e.g. Tacconi et al.|2006 ), only 
one of the SFRGs exceeds this limit, which is a dust opac- 



ity Eddington limit for SFR density (see Thompson et al. 
20051). This source is RGJ 123642, the highest redshift source 



in our sample having a low S/N in the radio, meaning the 
measured effective radius is unusually small. Even by assum- 
ing unresolved radio profiles, only RGJ 123642 exceeds the 
maximal starburst density, further highlighting that even 
starburst emission can dominate unresolved radio emission. 
In contrast, it would be unlikely for extended radio emission 
(this faint at these redshifts) to be driven by emission from 
AGN. 
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Table 3. Other Derived Properties of CO-observed SFRGs 
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Table Notes. Derived gas properties of SERGs. LpiR is derived from radio luminosity, and the associated SE R is thus labehed 
S'Fi^radio- Conversion to gas mass assumes the typical ULIRG conversion factor X =0.8 Mo (K kms~i pc^)"! (iDownes Solomonj 



||l998) for [U], the default assumption of gas mass for this sample, although we note that higher conversion factors 

(~4.5 Mo (K kms~i pc^)"!) might be appropriate for a subset of these galaxies which exhibit normal spiral galaxy characteristics (see 

[Daddi et al.|2010[ >, and these gas masses are given in [Sp]. The total dynamical mass is calculated by assuming a size of CO 

emission, which we do not directly measure, but assume is equal to the effective radius of radio emission, given as i^eff? cind when radio 

is not available we assume a radius of '^2.3kpc, the mean MERLIN+VLA radius, consistent with SMG sizes. The star formation 
efficiency (SEE) is S'Fi^radio divided by the ULIRG molecular gas mass, Mh2 • SER^/y has not been corrected for extinction. The star 
formation rate density, XIsfr, has units of Mo yr""*^ kpc"-^ {]). Since these galaxies were selected in the radio, we provide a measure of 
their starburst and AGN content by classifying them into starburst (SB), AGN, and starburst/AGN mix (SB+) in five separate AGN 

selection criteria. Class^^ marks any galaxy with X-Ray luminosity >5x 10^^ erg s"^ as "AGN", <lx 10^^ ergs"^ as "SB" and 
intermediate luminosit y gala xies as "SB+". X-ray data is taken from [Alexander et al.| ( [2003| >, |Brunner et al] ( |2008t , |Ueda et al.| ( [2008| >, 



and Mushotzky et al. ( 2000| >. Class^^ marks galaxies with significant 8/im fiux excess (relative to the stellar population fits) as "AGN", 
galaxies with marginally significant <2cr 8/im excess as "SB+" and galaxies with no 8/im excess as "SB" (see EigurejeJ. Class^^ 

marks 

galaxies with unusually high EIR-to-CO ratios as "AGN" {log{LFiR/ L'^q) >2.7), "SB+" (2.5 < log{LFiR/ L'^jq) <2.7), and at low 
EIR-to-CO ratios, "SB" {log{LpiR/ L'^q) >2.5). A fourth class may be defined as a measure of the star formation rate density, Ssfr; 
all of the galaxies for which we have Ssfr measurements satisfy the "SB" criteria however, with star formation densities 
<200 Mo yr""^- A fifth class is defined by their rest-UV/optical spectroscopic properties, which, by method of their selection, is 

uniformly "SB." 



The star formation efficiency (SEE) can be calcu- 
lated by dividing the FIR luminosity by the H2 gas mass. 
This calculation is contingent on the gas reservoir be- 
ing in the same region as the starburst, which is an as- 
sumption which needs to be investigated in more detail 
through future high-resolution multiple- J CO maps and re- 
solved FIR emission maps. The mean SEE for SFRGs is 
280di260 Lo/Mo- Put another way, the median depletion 
timescale (defined as Mhs/SFR, assuming 100% efficiency) 
for SFRGs is -^34Myr with inner quartile 20Myr - 55Myr. 
It should be noted here that a different gas conversion factor 
(e.g. 4.5 Mo (Kkms-^pc^)-^) would increase the depletion 
timescale and decrease the star formation efficiency by 5.6 
times. 



3.6 Dust Temperature and Dust Mass 

We fit FIR spectral energy distributions (SEDs) to the FIR 
fiux density limits at TO/xm, 350/im, 850/im, and 1200/im 
(see Table [l] for details). Our fitting method follows the 



methodology of Chapman et al. (2004), |Chapman et al. 
(2005|, and Casey et al. (2009a): a single dust temperature 
modified black body fit. We fix the dust emissivity (3 to 1.5 
and use the FIR/radio correlation to infer a FIR luminosity 
from the radio fiux. Since we are limited by a lack of data in 
the FIR for most of our sample (many only have limits), our 
derived dust temperatures are not well constrained. How- 
ever if the uncertainty in our assumption of the radio/FIR 
correlation and dust emissivity are taken as givens, then we 
deduce temperature uncertainties on the order of 20-30K. 
This is calculated by assuming Lfir scales with Lradi©; a 
higher Lradio implies a high Lfir, and a FIR SED with a 
high Lfir can then be constrained by the FIR fiux den- 
sity measurements, even in the case of upper limits. Despite 
the large uncertainties, our measurements (whose statistical 
mean is 65 K) are consistent with a selection of warmer dust 
galaxies, even though there may be outliers. 



We also estimate dust mass for the SERG sample by 
using the following relation 
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Mdust(M0) 



'S'ob 



(5) 



where Ki, is the dust mass absorption coefficient as 
a function of rest wa veiength: /i:850 = 0-15m ^ kg~^, 



K250 = 0.29=b0.03 m^kg-' jWiebe et ai.| |20Q9| ) and 
K7o = 1.2m^kg~^ jWeingar tner Dra ine||2001|). Again, the 
iack of FIR data means that our dust mass calcuiations are 
not weii constrained; our SFRGs have a mean dust mass 
upper iimit (at 2a) of (Mdust) <2xl0^ M©, which is not 
a weli constrained measurement. A quick comparison to 
the mean gas mass (2x10^° M©) impiies mean dust-to-gas 
ratios of <1/10 (weli within expectation). 



this quantity could be overestimated if is overestimated, 
the formation timescales likely represent a lower limit on 
the time it took to build up the stellar population since the 
star formation rates are hypothetically near their peak dur- 
ing the ULIRG phase ( ^ 2OOM0yr~^) in comparison to 
most galaxies at the same epoch (1-10 Mo yr"^). The hy- 
pothesis of SFRGs being near the peak is based on the fact 
that their SFRs are not sustainable beyond ^100 Myr, but 
it is unlikely that these systems are that young; it is more 
probable that SFRGs evolved more slowly, building up stel- 
lar mass gradually until a point, when a trigger led to an 
extreme star-bursting phase. 



3.7 Stellar Mass 

To estimate the galaxies' stellar masses, we combine the pho- 
tometric points in the optical (HST ACS B, V, i, and z 
for GOODS-N, HST ACS i for Lockman Hole, and Sub- 
aru i for SSA13) and the mid-IR {Spitzer-IRAC 3.6/im, 
4.5/im, 5.8/im and S.O/xm). For z =1-3, these photometric 
points cover rest-frame 1.6/im where stellar emission peaks. 
We use the hyperz photometric redshift code ( ^Bolzonella^ 



et al. 2000) to fit this photometry to several stellar pop- 
ulation SEDs fBruzual &: Chariot] [2003t. Using the best 



fit stellar population SED, we estimate the rest-frame H- 
band magnitude. While rest-frame K-band is often used to 
measure stellar mass for high-redshift galaxies (using the 
method outlined by Borys et al. 2005), recent work has 



shown that stellar masses derived from K-band are overes- 
timated due to increased contribution from AGN power law 
and decreased contribution from star light at longer wave- 
lengths, (see [Hainline et a"L]|2009 ). The mean light to mass 
ratio, Lh/M, can range from 5-10 Lq/M©, and using the 
'Bruzual Chariot] ( [2003| models, is set at 5.6 Lo/Mofor 
SMGs as in Hainline et al.. The mean rest-frame H-band ab- 
solute magnitude is -25.4zb0.8 with an inferred average stel- 
lar mass of 7x10^° M0. These masses are consistent with 
the H-band derived SMG stellar masses (H-band absolute 
magnitude averages -25.8±0.2 and masses average 2x10"^"^) 
found by Hainline et al. (2009 ) and a factor of ^2 lower than 
the K-band SMG masses derived by Borys et al. (2005). 



Figure [6] illustrates the best stellar population models 
with respect to the galaxies' photometric data. If an excess 
flux density is detected above the stellar model SED at ob- 
served 8/im then an AGN might be contributing significantly 
to the near infrared luminosities (see section ^ 



4.1 



for an anal- 
ysis of AGN contamination). While there are potential flux 
excess in seven of the 15 sources illustrated, the excess is 
only significant >2a in RGJ123653 a nd RGJ123718 . The 
former is discussed in detail in ^ Casey et al. ([2009a') who 
conclude an insignificant AGN contribution based on com- 
bined evidence from X-ray (very low flux, consistent with 
starbursts), extended radio emission (inconsistent with AGN 
radio emission) , and no clear powerlaw dominating the near 
to mid-IR data. The latter has extended radio emission, no 
X-ray detection and is the faintest radio galaxy of our sample 
(~15/iJy), and therefore, it is unlikely that it is dominated 
by a very powerful AGN. There is no clear relation between 
a system's stellar mass and its resulting detection in CO. 

The galaxies' formation timescales are given by Tform oc 
M^/SFR; the median Tform of the sample is lOMyr. While 



4 DISCUSSION 

Here we explore the relationship between the CO observa- 
tions in SFRGs and their other multi-wavelength proper- 
ties, drawing on comparisons with other galaxy populations. 
CO observations provide a unique and independent probe of 
the galaxies' star formation properties, by constraining the 
molecular gas reservoirs which fuel the star formation. In 
Table |4] we have summarised many of the measured phys- 
ical properties of SFRGs relative to SMGs for quick refer- 
ence. Important to note, however, is that the SMGs of Ta- 
ble [4] are those with CO observations, limited to the mostly 
bright SMG subsample. After comparing with other popu- 
lations, we discuss the implications that these observations 
have for the study of all high-z ULIRGs and how improved 
targeted observations, in both molecular interstellar medium 
lines and FIR continuum, will enable thorough, unhindered 
analysis of ULIRG evolution and extreme star formation at 
z-2. 

The original motivation for segregating submm-faint 
and submm-bright populations for CO observations is the 
premise that they exhibit similar extreme starburst qual- 
ities, thus molecular gas qualities, but differ in dust dis- 
tribution whereby slightly warmer dust systems are unde- 
tectable at 850/im (the usual SMG-selection band). Survey- 
ing SFRGs in CO provides confirmation, through detection 
of vast gas reservoirs, that SMGs are not the only significant 
population of ULIRGs at high-z. The basic physical premise 
of the comparison is that the dust distribution in SFRGs 
would need to be clumpier or more compact to be heated to 
slightly higher dust temperatures. Since SFRG radio mor- 
phologies on average seem to be similarly extended and ir- 
regular, SFRGs are suggested to be less homogeneously dif- 
fuse than S MGs ( [Menendez-Delmestre et al.||2009| |Hainline| 
et al.|2009t , with more concentrated clumps, yet spread over 
the same large area. 



4.1 Mult i- wavelength Properties of SFRGs and 
AGN fraction 

Table |3] provides a summary of SFRGs' physical prop- 
erties derived from multi-wavelength data. It includes 
AGN/starburst classifications for each individual SFRG, 
which help shed light on the complex nature of the popu- 
lation, and the potential for AGN contamination. "Class'^" 
classifies SFRGs based on X-ray luminosity, "Class^^" clas- 
sifies according to S/xm flux excess in the near-IR, and 
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Figure 6. Stellar population models are fit to SFRG optical and near-IR photometric points. We estimate stellar mass using the 
rest-frame H-band magnitudes which are inferred from the fits; the method is described in section [3?f| The SFRGs which are detected 
in CO have their names enclosed in boxes. 



"Class^^" classifies according to FIR-to-CO luminosity ra- 
tio (where an unusually high ratio can be accounted for 
by an AGN contributing significantly to radio luminosity, 
thus overestimated FIR luminosity) . A fourth class and fifth 
could also be defined, as rest-UV/optical spectroscopic fea- 
tures and as a function of the star formation rate density as 
traced by MERLIN+VLA morphologies, however we have 
noted that all of the samples' morphologies are consistent 
with starbursts in both criteria. While we leave the reader 
the flexibility to interpret the AGN content of the SFRG 
sample as (s)he sees flt, we infer a total SFRG AGN frac- 
tion of O.SdiO.l to first order (included in Table^, the mean 
AGN dominated fraction from each of the three Table [3] clas- 
sifications. This agrees with earlier measurements and esti- 



mates for SFRGs given in Casey et al. (2QQ9bJ and Casey 



et al. (2009a). Also, the possibility exists that weak beaming 



of low- luminosity radio jets could overestimate radio lumi- 



nosity by factors of a few for some SFRGs (please see Casey 
|et al.||2009b[ for the detailed discussion of radio beaming in 
more compact SFRGs). Again, we caution AGN contami- 
nation has the potential to strongly bias the CO observa- 
tions/interpretation of individual objects. 

Figure [5] shows that SFRGs have faint optical lumi- 
nosities, suggestive of heavy reddening or obscuration by 
dust. This idea is supported further by the large discrep- 
ancy between star formation rates derived from radio lu- 
minosity versus rest-UV flux densities. Even after correc- 
tion for dust extinction using the UV slope and dust ex- 
tinction models of |Calzetti et al. ( |1994| , the UV derived 
SFR can be a factor of 10-100 times lower than the radio- 
inferred SFR demonstrating that SFRGs are far dustier than 
normal dust extinction laws predict (which also holds for 
SMGs). Their radio morphologies are primarily extended 
and irregular, thus attributable to star formation and not 
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Table 4. SFRG Properties Summarised Relative to SMGs 



Property 


SFRGs 


SMGs" 


iz) 


1.8±0.7 


2.4±0.6 


IK V J 


(4±3)xl0-'^^ 


(8±4)xl0-'-^ 




700±500 


1400±500 




2.3d=0.8 


2.7d=0.4 




30+^0 

— 20 




ausTj \ / 


66±15 


41±5 


Mdust (Mo) ^ 


<2xl09 


9x10^ 


H-band Mag 


-25.4±0.8 


-25.8±0.2 


(Mo) 


(7±6)xl0io 


(2±l)xl0ii 


{rform) (Myr) 


lot- 




Mh2 (Mo) 


(2.1±0.7)xl0^'-' 


(5.1±2.8)xl0^° 


AVco (kms-i) 


320±80 


530±150 


SFE (Lo/Mo) 


280±260 


450±170 


Mdynsm^i (Mo) 


(7.2l^-j)xlOio 


(l.SlJDxlOii 


(Tdepl) (Myr) 


34±24 




fgas = (Mgas/Mdyn) ^ 


-0.02 


^•^^-0.07 


Utars = {M^/Mdyn) 


-0.2 


-0.3 


'S'cO(4-3)/'S'cO(3-2) 


2.0±0.8 


1.4±0.2 


/agn ^ 


0.3±0.1 


0.4±0.2 



" All aggregate properties of SMGs are measured with respect 
to the CO-observed subset from |Neri et al.| ( |2003| >, [Greve et al.| 
( |2005t , and |Tacconi et aLlpOOef . 

^ Effective radius and SFR density are measured for MER- 
LIN +VLA_^mage£^SFRG^ only. The effective size of SMGs is 
from |Biggs Ivison| ( [2008t and the SFR density is the SMG SFR 
divided by the mean Rgff • 

Tdust cind Mdust fits for SFRGs are described in section |3.6| 
while both measurements are highly uncertain due to poor FIR 
flux density constraints, their calculation is useful for comparison 
with SMGs from |Chapman et al.| ( |2005t and |Kovacs et al.| ( [2006t . 
^ Dynamical mass assumes an effective radius of 2 kpc for SFRGs 
and SMGs without radio size measurements; for SMGs, this size 
is supported by measurements from [Tacconi et al.| ( |2008[ ) , and 
for SFRGs the size has been shown to be similar for one SFRG 
( [Bothwell et al.|2010| and is similar to their MERLIN+VLA radio 
sizes (Reff). 

^ The gas and stars fraction represent the fraction of total mass 
which is in gas (or in stars), internal to each galaxy. Note: the gas 
fraction does not include the 40% correction for helium. 
f The ratio of CO line fluxes represents the source excitation. 
The value for SFRGs is based on the single measurement of 
RGJ123711 described in section [3^ and the SMG measurement 
is taken from the three SMGs measured by (Weiss et al.| ( |2007^ ). 
^ The AGN fraction (within a population) is estimated in the 
SFRG sample as described in section [4.31 while the AGN fraction 
of CO-observed SMGs is taken from their rest-UV spectral classi- 
fication (note none of the SFRGs have AGN spectral signatures). 
The entire SMG population is estimated to have an AGN fraction 
of -0.25 (see (Alexander et al.|2005t . 



compact AGN as is also the case for SMGs fBiggs k Ivi- 
son 2008) [Chapman et al.|[2004|). We also note the effec- 



tive MERLIN+VLA radii between SFRGs that have been 
CO detected (i?eff = 2.7±0.6 kpc) and those that have not 
{ReS = 1.7±0.8kpc). While the difference is insignificant, it 
could be explained by observing that the sources with lower 
effective radii have lower S/N in the radio, implying lower 
FIR luminosities. If we assume the inferred radio/FIR and 
the ULIRG L^q/Lfik relations, we predict integrated CO 
fluxes of 0.1-0.3 Jy kms~^ for CO-undetected sources based 



on radio flux density (requiring particularly sensitive obser- 
vations with our setup). 

Figure [6] showed that most SFRGs have only minor 
AGN contribution in the near-IR; most emission is dom- 
inated by stars. The stellar masses of SFRGs aver age to 
7x10^° Mo, consistent with SMG masses derived by 



Main- 



line et al. ( 2009 ) using rest H-band luminosity but a factor 



of ~2x lower than the SMG masses of Borys et al. (2005) 
from rest K-band luminosity. 

4.2 Comparing to SMGs 

The comparison of molecular gas properties between SFRGs 
and SMGs is the main focus of this work, however, a 
strong selection bias is folded into this comparison. The 



first SMGs to be observed in CO (Neri et al. 



2003 



Greve 



et al. 2005 Tacconi et al. 2006) were among the bright- 



est SMGs with spectroscopic redshifts (Lfir ^10 " Lo) 
since brighter Lfir systems had a higher likelihood of be- 
ing CO detected. We note that the luminosity distribution 
of SMGs plotted in Fig. [l] includes many lower luminosity 
SMGs which have recently been observed in CO, in par- 
allel to our SFRG observing programs, but have not yet 
been analyzed or published (Bothwell etal., in preparation). 
Most of the 'bright-SMGs' (those with published CO spec- 
tra) have AGN signatures in their rest-UV spectra, a prop- 
erty which would exclude them from SFRG selection if they 
were submm- faint. This means that similarly bright radio 
sources which are submm-faint, comparable to the 'bright- 
SMGs,' were excluded from our sample due to AGN con- 
tamination. 

Besides the luminosity bias introduced by weeding out 
AGN as revealed by rest-UV/optical spectroscopy, a fur- 
ther bias exists due to spectroscopic incompleteness of the 
SFRG population. Since their radio emission was more likely 
thought to be dominated by AGN, SFRGs were not fol- 
lowed up in rest-UV/rest-optical spectroscopy nearly as 
thoroughly or completely as SMGs. This likely means that 
the absolute brightest SMGs were CO observed while a large 
sample of bright SFRGs could have been excluded from CO 
observations due to a lack of reliable redshift information 
or potentially strong AGN contamination. Furthermore, the 



original SFRG selection of Chapman et al. (2004) had the 
added criteria of a faint optical magnitude i > 23, which 
made redshift measurement from rest-UV spectra more dif- 
ficult. There has been some anecdotal indication, however, 
that warmer dust systems might not exist at the highest lu- 
minosities with the same frequency as colder dust systems, 
as revealed by 250/im selected HyLIRG populations ( Casey] 
eraL|[20TQ| . 

We note that the mis-identification of optical counter- 
parts to SMGs has potentially led to a lower CO-detection 
rate for that population than other CO-observed galaxy pop- 
ulations; some of the radio galaxy SMG-counterparts might 
not be starbursts and might have intrinsically low FIR lu- 
minosities, thus CO luminosities. We note that SFRGs (al- 
though not faced with the issue of matching FIR positions 
to a correct radio counterpart) could also suffer from the 
selection of non-starburst radio galaxies. 

In the sections below, we frequently discuss how SFRGs 
relate to SMGs: both the CO-observed 'bright' subsample of 
SMGs and extrapolations based on preliminary analysis on 
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the fainter, more numerous sample of SMGs (Bothwell, pri- 
vate communication). The overall luminosity bias existing 
in these distinct samples must be kept in mind when popu- 
lation comparisons are made. 

4.3 Star Formation Efficiencies 

The relationship between star formation rate and molecu- 
lar gas mass is paramount to a galaxy's evolutionary in- 
terpretation. This is measured by comparing the FIR lu- 
minosity with CO line luminosity, as we show in Figure jT] 
In this context, SFRGs appear to have similarly high CO 
luminosities as SMGs, and most SFRGs lie slightly above 
the 'ULIRG' star formation efficiency powerlaw relation, 
Leo Lfir°"^^, which is followed by both local ULIRGs 
and SMGs. However, SFRGs are inconsistent with the star 



galaxies {L'qq oc Lfir°"^^; Solomon & Sage 


1988 Gao & 


Solomon 2004). The SMGs shown on Fig 


7 


are those of 


Neri et al. (2003), Greve et al.| (|2005| and Tacconi et al. 



(2006}. As mentioned in section |4.2[ SFRGs are a factor of 
~2 less luminous in radio (thus Lfir) than CO-observed 
SMGs. Due to their relative high luminosities with respect 
to local ULIRGs, the SMGs have been described as the 



scaled-up high-redshift analogues to local ULIRGs (Tac- 
coni et al. 2006 ). In contrast, SFRGs probe a luminosity 
regime ^10 ' Lq closer to the locus of local ULIRGs at 
~10^^ L0. While SFRGs might seem to be better analogues 
of local ULIRGs than SMGs in luminosity space, we note 
that the stellar and gas properties of the two populations 
are quite distinct: local ULIRGs being more compact with 



lower stellar masses than SFRGs and SMGs (Dasyra et al 



2006), and they have lower CO luminosities by a factor of 
-2-3. 

The star formation efficiencies (SFEs) of SMGs and 
SFRGs span the range 70-1000 L©/ M© . The median SFE of 
the CO-detected SFRG sample is 280±260 L© M©"^ which 
is statistically consistent with the mean SFE for SMGs, 
450±170 L© M©~^, although both values incorporate the 
large uncertainties of the gas conversion factor and FIR- 



derived SFR. Chapman et al.^ ( ^2008 J highlights that the two 
pilot program '^CO(J=3^2) detections of RGJ123711 and 
RG J 163655 have exceptionally high SFEs and hypothesised 
that SFRGs, with further observation, might show similarly 
high SFEs compared to SMGs. |Daddi et al.| ( 2008 ) analyzed 
the CO content of two BzK selected galaxies (RG J 123710 
and RGJ123626, also selected as SFRGs and included in our 
analysis) and claimed that they had relatively low. Milky 
Way/ "normal spiral" efficiencies, emphasising the difference 
between them and the high-efficiency ULIRGs. Our large 
sample of SFRGs, including both the Chapman et al. and 
Daddi et al. subsamples reveals a much wider spread in star 
formation efficiencies, suggestive of a wide range in gas states 
and a possible range of galaxy states, although more SFRGs 
are consistent with the less- efficient Daddi et al. sample. 

The SFRGs not detected in CO and those far below the 
ULIRG Lpm/L^Q relation would appear to be very efficient 
star formers (less gas to fuel their high SFRs), however this 
assumes that AGN contamination is minimal. AGN contam- 
ination is more likely than super-efficient star formation and 
happens when AGN boosts the radio- inferred FIR luminos- 
ity, thus inferred star formation rate; as the empirical rela- 
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Figure 8. The distribution in CO line widths for the SFRG sam- 
ple compared to local ULIRGs ( [Solomon et al.|l'997j and SMGs. 
The offset/tentative SFRGs (open histogram) are added on top 
of the the remaining SFRG sample (hashed area). The distribu- 
tion for local ULIRGs has been re-normalised with respect to the 
total number of SFRGs for a more clear comparison. The SMG 
distribution has been corrected for overestimation in line widths 
which is caused by fitting a single Gaussian to a double-peaked 
CO line, and the SMG distribution includes the samples of Neri] 
leTal-l ( |2003| , [Greve et al. (2005), and Tacco ni et al.[ |2006), as 
well as some yet unpublished SMG observations (Bothwell etal., 
in preparation). 



tion between Lfir and Lqq suggests, AGN contaminated 
sources would be fainter in CO gas than predicted. We then 
can classify SFRGs in terms of CO luminosity to FIR lu- 
minosity (i.e. the ratio of Lco/^fir), where low ratios are 
designated 'AGN' in 'CLASSco' in Table|3] While the scat- 
ter of local ULIRGs around the Lqq/Lfik relation is mini- 
mal, '-^O.Sdex, both SFRG and SMG populations have more 
significant scatter below the relation; this is consistent with 
many SFRGs and SMGs having powerful AGN which boosts 
radio (thus FIR) luminosity. 

The scatter of SFRGs above the ULIRG Lco/^fir re- 
lation is suggestive of low star formation efficiencies as pre- 
sented by [Daddi et ah] ([20081 and [Daddi et al"! ([20101 . We 



caution that sources like RGJ123710 and RGJ123626 which 
have these low star formation efficiencies, it might be more 
appropriate to assume that the gas properties are more sim- 
ilar to spiral/disk galaxies than ULIRGs, particularly when 
converting to molecular gas mass using the CO/H2 conver- 
sion factor. This factor differs greatly between ULIRGs (0.8 
Mq I (K kms~^ pc^)) and Milky Way type disk galaxies (4.5 
M0/(K kms-^pc^)). Since the inferred gas masses differ so 
greatly given these different assumptions, we give for 
both ULIRG and spiral/disk galaxy assumptions in Table [s] 
but proceed with our interpretation using the ULIRG con- 
version factor. 



4.4 Line Widths: Implications for Merger Stage 

The distribution in CO line widths provides important in- 
sight into the galaxies' dynamics. Figure [S] shows the AVco 
full width at half maximum (FWHM) distributions for our 
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Figure 7. Left: FIR Luminosity against CO Luminosity, ^co\i-o ]' ^^^^^ {black c ircles) lie in the same luminosity space as local 
ULIRGs {gray trmngr/eg; [Solomon et al.|1997| >, while SMGs (croggeg; [Greve et al. [20051 have higher luminosities. SFRGs detected in CO 
are solid while SFRGs without CO detection are open circles. The tentative identifications or "oflPset" CO sources are circles with gray 
centres. We also overplot data of spiral galaxies for perspective {gray dmmon(ig; [Solomon Sage|1988|[Gao Solomon|2004| ), and include 
the best-fit observed relations between LpiR and LJ^q for local spirals (dotted line) and local ULIRGs/SMGs (dashed line). Right: the 
star formation efficiency, given by Lfir/Lco? is plotted with redshift. SFRGs appear to share the same range of SFEs as SMGs despite 
their fainter luminosities. The mean SFE and 1-a bounds of local ULIRGs are illustrated by the horizontal solid and dashed lines. 



SFRGs, SMGs (both 'bright' SMGs and the fainter sample 
observed in CO only recently; Small, private communica- 
tion) and local ULIRGs (Solomon et aL[[1997|). None of the 
line widths presented here have made any inclination an- 
gle assumptions. The line widths of SMGs were adjusted to 
correct for the prior exclusion of double-peaked Gaussians 



(see the details of this correction in Coppin et al. 2008); 
this has reduced the mean SMG line width from 600 kms~^ 
to 530 kms~^. Despite the adjustment, SMGs seem quite 
distinct from SFRGs and local ULIRGs by having a high- 
FWHM tail in its distribution. This high-FWHM tail is seen 
only in bright subsample of SMGs originally surveyed in 
CO gas. While it could be attributed to selection bias, in 
that wide CO features are only detectable in the brightest 
objects where S/N is much higher, we highlight that there 
was significant improvement in receiver sensitivity between 
observations of these bright and wide SMG CO lines and 
the fainter SFRG observations, so the data have compara- 
ble S/N. This raises the possibility that only the brightest 
subsample of high-z ULIRGs (Lfir ^ 10^"^ Lq) have wide 
CO line widths ( > 500 kms"^). 

Broad dispersion-dominated CO lines (and irregular 
double peaked profiles) in the highest luminosity systems is 
suggestive of early stage major mergers where two gas rich 
disks are infalling. Local ULIRGs in contrast have a much 
narrower line width distribution and are a factor of ^-^5-10 
fainter in Lfir. For this reason, local ULIRGs are often said 
to be in a late starburst phase, at a coalesced point during 
a merger (when progenitors have coalesced into a single sys- 
tem; for a review see [Sanders k Mirabel] 1996| ). SFRGs and 
more modest-luminosity SMGs are difficult to place in this 
evolutionary sequence, but their populations are not likely 
to be exclusively dominated by either beginning or ending 
merger sub-stages. 



The lower dynamical masses of local ULIRGs (typical 
sizes R = 1 kpc and AV =300 kms~^) could be due to 
downsizing— where extreme starbursts today are less mas- 
sive than those at high- 2;. Both SFRGs and SMGs are con- 
sistent with this picture since they seem to be a factor of 
^2 larger {R1/2 ^ 2 kpc, see measurements of CO size in 



[Tacconi et al.|2008[[Daddi et al.|2010[ |Bothwell et al.|2010 ) 

The subsample of SMGs which have very broad features may 
be in a particular stage of merger where their line profile 
becomes broadened, perhaps in observation of two distinct 
components with very small physical separation or gas. This 
would appear to increase our dynamical mass estimates of 
these systems. However, the radial 'size' and merger cor- 
rection factor C would both need to be reasessed (both of 
which have not been measured on a case by case basis for 
these sources) before physical interpretation of higher dy- 
namical mass estimates were made. 

Although more observations need to be taken to con- 
clude, an alternate explanation for the narrower line widths 
observed in modest-luminosity sources, including SFRGs, is 
that the population consists of fewer major mergers than the 
very bright systems. Observational evidence indicates that 
anywhere from 50-90% of local ULIRGs have undergone re- 
cent mergers (Lawrence et al.|1989 Melnick & Mirabel 1990 



Clements et al. |1996| e.g.), but that a sizable fraction might 
be triggered by other mechanisms. Recent work from Gen- 



zel et al. (2008) suggests that high star formation rates in 



secularly evolving disk galaxies may be caused by rapid ro- 
tation or smooth accretion of material from its surroundings 
(e.g. minor merging or tidal accretion). In addition, theoret- 
ical work indicates that the high star formation rates and 
IR luminosities in ULIRGs could often be generated by mi- 
nor mergers and turbulent disk processes (e.g. Monaco|2004 



Dib et aLf2006J . ULIRGs driven by secular processes would 
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exhibit narrow CO line widths, consistent with the SFRGs 
presented in this paper. High spatial resolution gas obser- 
vations are needed to determine the true nature of their 
dynamics however. 

Recent observations and simulation work have reiter- 
ated the idea that >10^^ ULIRGs are more often in early- 



stage major mergers. Tacconi et al. (2008) and EngeletaL 
( |2010| show that most SMGs at z~2 exhibit disturbed gas 
morphologies rather than smoothly rotating disks, and sim- 
ulations and semi- analytic SMG models tell us that major 
mergers likely initiate most ultraluminous phases of high-z 
star formation seen in SMGs ([Narayanan et al.|2009 



bank et al.||2008 



from Dave et al 



Baugh et al. 2005| . However, recent work 



(2010) indicate that ULIRGs may also be 



driven by continual bombardment by very low mass frag- 
ments onto a '-^lO^^ M0 galaxy. Cold streams feeding con- 
tinual gas buildup (e.g. Dekel et al. 20091) has also been 



raised as a possible origin. The tail of large CO line widths 
in SMGs provides a crucial piece of observational evidence 
that some SMGs are much more highly disturbed and rep- 
resent a different phase than SFRGs and local ULIRGs. 

We note that the mean gas fraction, defined as 
gas mass to dynamical mass ratio, of SFRGs is 
(Mgas/Mdyn)= /gas= 0.071q J2i cousistcut with the same ra- 
tio for SMGs, which have (Mgas/Mdyn) -O.OQIq:?? (after 
correction for a 30^ inclination angle). While this compar- 
ison is between SFRGs which have lower luminosities than 
CO-observed SMGs, the same molecular gas fraction sug- 
gests that the two populations are likely in similar evolution- 
ary stages. We note, however, that if we assume a CO/H2 gas 
conversion factor consistent with spirals instead of ULIRGs, 
the gas mass fraction increases substantially to ^0.60. While 
overall, SFRGs and SMG properties are more consistent 
with ULIRGs, it is possible that a few outliers, for exam- 
ple RGJ123710 and RGJ123626 described in fPaddi et al 



2008 ) and a few of the SMGs exhibiting unusually low SFEs, 



are much more gas rich than the ULIRGs which comprise 
the rest of the populations. 

4.5 Comparison with Other Populations 

An evolutionary transition stage between ULIRG and 
quasar has been explored extensively by work on Dust Ob- 
scured Galaxies (DOGs; Dey et al. 20081 [Pope et al.|| 2008 |), 
a population of galaxies with warm dust temperatures and 
more modest SFRs than the brightest SMGs, consistent with 
-10^^ L© ULIRGs. Ten of fourteen SFRGs with 24/im flux 
measurements satisfy DOG selection. While DOG selection 
is quite broad and selects nearly all SFRGs and many SMGs 
(requiring 24/im flux densities >100/iJy and red optical to 
IR colors), a subset of DOGs with spectroscopic redshifts 



have detailed near-IR to FIR photometric constraints ( |Buss- 
'mann e t al.|2009 ) which show that AGN may contribute sig- 



nificantly to their bolometric luminosities, and as a result, 
most have warm dust temperatures {Td ^ 45 K). Due to 
their high AGN fraction (dependent on luminosity and only 
^0.5 in the faintest, S24 <0.5mJy subset, e.g. [Pope et al.| 
[200 8 ) and heavy dust obscuration, many DOGs are believed 
to lie at the transition phase between SMGs (or star-forming 
ULIRG) and luminous quasar ( |Pope et al.|2008 ) and overlap 
with the SFRG population. 

While SFRGs might have warm dust temperatures like 



10-'' 




? SFRGs 
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Figure 9. The star formation rate per unit stellar mass 
against stellar mass. We compare SFRGs {large circles) to SMGs 
(crosses) and z~l starburst galaxies (gray triangles and squares; 
[Noeske et al.ll^OOTl ). We also overplot the derived redshift de- 
pendent relations (at z=l and z=2, dot-dashed lines) found for 
GOODS galaxies in paddi et al.| ( |2007| >. SFRGs not detected in 
CO are open circles while detected SFRGs are filled. Like SMGs, 
SFRGs have very large star formation rates per stellar mass com- 
pared with "blue sequence" galaxies, although their mean stellar 
mass (~7xlOio M0) is less than the mean SMG stellar mass 
(~2xl0ii M0). 



some DOGs, we find that most are dominated by star forma- 
tion and not AGN. This is largely a function of the aggres- 
sive selection of SFRGs, meant to weed out strong AGN by 
their spectral indicators in the rest-UV/optical and the ob- 
servation of minimal 8/im flux excess and of extended radio 
emission. High star formation rates and a low AGN frac- 
tion (with respect to a higher AGN fraction in DOGs) are 
strong evidence that SFRGs are at a similar ULIRG phase 
to SMGs despite their warm dust. 

The molecular gas fractions of normal starburst 
galaxies have been estimated at /ga s ^^0.4-0.5 (assuming 



^co/ L 



FIR prior; 



Erb et al. 2006). The gas fraction in 



SMGs (and a few BX active galaxies) has been measured 
to be /gas ~0.3-0.5 (these values do not take inclination 
into account, which is a factor of 1/4 [Tacconi et al.||2^08| . 
We measure an internal gas fraction of our SFRG sample 
of ^0.07, which is consistent with SMGs, both the Tac- 
coni et al. estimate (/gas ^0.1 if corrected for inclination) 
and our reassessment of the same data (/gas =0.091° °^). 
All /gas measurements for high-z ULIRGs (/gas ^^0.1) are 
lower than the inferred gas fractions for normal z^2 galax- 
ies (/gas ^^0.5) from Erb et al. (2006). This could indicate 
that the more modest-luminosity galaxies, consisting of gas- 
rich disks, have low dynamical masses but proportionately 
high gas mass, meaning they would be good progenitor can- 
didates for ULIRG systems, if set on collision courses with 
other gas-rich disks. During the ULIRG starburst phase, the 
gas mass would start to decrease with rapid star formation. 
In contrast, we recognize that by using a higher Xco con- 
version factor more consistent with quiescent disks, then the 
measured gas fractions of these ULIRGs would increase by 
'-^6x, making their gas fractions consistent with those of 
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the estimate for normal z ^2 galaxies. Future observations 
of more modest-luminosity galaxies in CO gas (e.g. like the 
recent work of Tacconi et al.|l2010) are needed to to truly 
understand the sequencing and gas properties of ULIRGs 
and their progenitors. 

Figure [9] highlights the unusually high star formation 
rates per unit mass of SFRGs and SMGs above normal star- 
bursts (AEGIS samples; [Noeske et al.|2007 |), for a wide range 
of stellar masses. The relation between stellar mass and spe- 
cific star formation rate has been shown to evolve with red- 



shift, e.g. Daddi et al. (2007); however, SMGs and SFRGs 



still lie at higher star formation rates than galaxies of equal 
mass at the same redshift (e.g. |da Cunha et al. 2010). This 
enhanced SFR fraction per unit mass reveals that SFRG and 
SMG star formation processes seem fundamentally different 
than SF processes in more modest luminosity galaxies, de- 
spite the overall range of stellar mass spanning almost two 
orders of magnitude. 

The BzK active galaxy selection ( [Daddi et al.|2004| is 
meant to select moderate star forming (~1OO-2OOM0 yr~^), 
massive (^^10^^ M©) galaxies at high redshift— systems 
which are typically below t he ULIRG star for mation rate 
threshold ( > 200 Mq yr~^).|Daddi et al.|(|2010t detect sev- 



eral active BzK galaxies in C0( J=2^1) with surprisingly 
high gas masses given their star formation rates, indicating 
that they exhibit gas properties more consistent with nor- 
mal Milky Way type galaxies, but at 'scaled-up' luminosi- 
ties, stellar masses and gas masses. It is important to note 
however that the active BzK galaxies observed in CO were 
all radio-detected, in other words, most of them have SFRs 
above the ULIRG cutoff (>2OOM0yr~^) and might other- 
wise be characterized as SMGs or SFRGs. All SFRGs in this 
paper satisfy the active BzK selection criterion, indicating 
that BzK selection might probe both massive star-bursting 
galaxies and massive, extreme, dusty starbursts. Similarly, 
it appears as if SFRG selection might select both high-z 
ULIRG merging systems and extreme gas-rich disk galaxies. 



4.6 Volume Density 

It is probable that the luminosity bias described in sec- 
tion |4.2| has significant effect on how we can interpret the 
CO observations of either population. 

Spectroscopic incompleteness in the SFRG population 



(as discussed in section 4.2) makes volume density diffi- 
cult to calculate. We estimate that about 35 % of a com- 
plete sample of /iJy radio galaxies (^0.7arcmin~^) with 
Si. 4 >20/iJy have not been observed spectroscopically (see 
Chapman et al.|2003 ). The majority of these are submillime- 
tre faint (since most spectroscopic observations of //Jy radio 
galaxies has been of SMGs). Of the submm- faint galaxies 
which were spectroscopically observed, 65 % have confirmed 
redshifts, roughly half of which are low luminosity AGN and 
the other half are star forming galaxies. This means that 
by completing spectroscopic observations of known submm- 
faint /xJy radio galaxies in well surveyed fields, the source 
density of star-forming, submm-faint ULIRGs with redshifts 
would likely increase by ^-^30 %, increasing the source density 
of all high-z ULIRGs by ^ 15%, which could also increase 
the ULIRG contribution (including SMGs) to the cosmic 
star formation rate density at its peak (e.g. see 
et al.||2009||Goto et al.||2010| ). 



Bouwens 



Future observations from SCUBA2 and the Herschel 
Space Observatory in the 50-500 /im range will dramatically 
improve the census of warmer dust ULIRGs (and all dusty 
starburst galaxies) at z^2 through continued discovery, thus 
allowing more thorough follow up of the submm-faint (and 
50-500/im bright) radio sample. With improved interfero- 
metric millimetre line observations at low and high frequen- 
cies from the Atacama Large Millimeter Array (ALMA), we 
will be able to target these high redshift sources in multiple 
CO transitions, thus removing an observational bias towards 
certain J-transitions of CO and enabling more accurate cal- 
culation of gas masses. Without such strong temperature 
biases in gas and dust observations, a more complete inter- 
pretation of high redshift ultraluminous galaxies will finally 
be possible. 



5 CONCLUSIONS 

We have presented CO molecular gas observations of a 
sample of submillimetre-faint, star-forming radio galaxies 
(SFRGs). Due to their non-detection at submillimetre wave- 
lengths and lack of dominant AGN, these ultraluminous, /xJy 
radio galaxies are thought to be dominated by star forma- 
tion but have warmer dust temperatures than SMGs. Out 
of 16 CO-observed SFRGs (12 from this paper and 4 from 
the literature), 10 are detected with a mean CO luminosity 
of 

^co[i-o] 2.6x10^° K km s ^pc^, which is slightly less 
luminous than the CO-observed SMG sample, despite being 
r^2x less luminous in the radio. We attribute the luminos- 
ity difference to a selection bias but suggest that physical 
driving mechanisms might differ between the very bright 
(>10^^ L0) and moderately bright (^10^^ L©) populations. 

High-resolution radio imaging from MERLIN+VLA 
shows that the radio emission in the SFRG sample is re- 
solved and extended with mean effective radii ^2kpc, sug- 
gesting that the SFRG radio luminosities are dominated 
by star formation rather than AGN. The MERLIN+VLA 
sizes constraints are consistent with similarly analyzed SMG 
MERLIN+VLA sizes. While we note that AGN do not dom- 
inate our sample, it is possible that several of our sources 
have non-negligible AGN due in part to their selection as ra- 
dio galaxies. Due to limited FIR data, we use the FIR/radio 
correlation to derive L fir and then compute extinction- free 
star formation rates from the FIR. 

The star formation efficiencies (SFEs) of SFRGs are 
comparable within large uncertainties to those of SMGs and 
local ULIRGs, even though a few sources appear to have 
very high SFEs (like tho se in Chapm an et al. 2008 ) or very 
low SFEs (like those in jDaddi et al.|2008t . Those with per- 
ceived very high SFEs are more likely AGN dominated than 
super efficient; their FIR luminosities as calculated from the 
radio are probably overestimated. 

SFRGs have narrower CO line widths than the bright 
subsample of SMGs at the same redshifts (AVsfrg ^ 
320 kms~^ and AVsmg '^530 kms~^), suggesting that 
SFRGs might have less disturbed dynamical environments. 
The line width distribution is potentially suggestive of dif- 
ferent evolutionary stages or processes between SFRGs and 
SMGs; however, the observed difference with SMGs could 
be due to a S/N or luminosity bias. The former would mean 
that intrinsically broad lines would have underestimated 
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FWHMs due to low S/N. The latter is due to the more thor- 
ough spectroscopic sampling of the SMG population. SMGs 
have higher spectroscopic completeness and also include 
many objects with AGN signatures in the rest-UV/optical. 
Any SFRGs which have similar AGN signatures were culled 
from the sample, thus eliminating some of the potentially 
brightest SFRGs (most of the bright SMGs which have been 
surveyed in CO contain optical AGN). While less luminous 
SMGs (at the same luminosities of SFRGs) exist, few have 
been observed in CO, thus it is difficult to rule out that CO 
luminosity might relate to the distribution in CO with line 
width. 

Despite selection biases, we have explored the possible 
physical scenarios triggering warm-dust ULIRGs in contrast 
to the well studied cold-dust ULIRGs. SFRGs appear to 
bridge the gap between the properties of >10^'^ LqSMGs 
and ~10^^ L0. Quantitatively, their extended radio emis- 
sions suggest sizes consistent with SMGs, implying much 
larger dynamical masses than local ULIRGs. We show that 
SFRGs have the same AGN fraction as SMGs and are there- 
fore unlikely to represent a 'post-SMG' AGN turn-on phase. 
Luminous SMGs have been characterised as an early infall 
stage during a major merger, and local ULIRGs are often 
described at late- type major mergers. Here, we suggest that 
SFRGs (and less luminous SMGs) span the range of states 
during peak merger interaction. 
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